Measurement of helium-3 and deuterium stopping power ratio for negative muons 
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The measurement method and results measuring of the stopping power ratio of helium-3 and 
deuterium atoms for muons slowed down in the D/^He mixture are presented. Measurements were 
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experiment was carried out at PSI muon beam pE4 with the momentum P/i = 34.0 MeV/c. The 
measured value of the mean stopping ratio Sajfe/D is 1.66 ±0.04. This value can also be interpreted 
as the value of mean reduced ratio of probabilities for muon capture by helium-3 and deuterium 
atoms. 
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I. INTRODUCTION 

Atomic capture of pions and muons stopped in a 
mixture of hydrogen and helium isotopes has been 
a subject matter of quite a lot of experimental Q 
12, B , B , p. B , l a B and theoretical investigations 
llfllll 111 ElIIlll III III El El 13. The inves- 
tigation of this process is important for understanding 
physics of exotic system and also for studying elemen- 
tary processes occurring when negatively charged par- 
ticles stop in a material. To separate two processes - 
atomic capture and transfer of muons from /i-atoms of 
hydrogen isotopes in the course of their de-excitation 
to helium nuclei - is practically impossible if x-rays 
are used, as usually, as diagnostics. What one usual 
observes experimentally is the result of interference of 
a few processes accompanying the muon capture by 
atoms of hydrogen and helium isotopes. Therefore, it 
is quite a problem to extract unambiguous informa- 
tion on the law of initial capture of muons by atoms 
of the hydrogen (deuterium)-helium mixture compo- 
nents. The earlier made assumption that the prob- 
ability for direct muon capture by atoms of the me- 
chanical H2(D2)/M (M is 3^4jjg^ Aj.^ mixture is 
proportional to the charge and concentration of each 
component (Fermi- Teller Z-law turned out to be 
wrong. Experimental data E3 revealed deviation from 
the Fermi-Teller law. The Z-law is based on the as- 
sumption that the atomic capture probability is pro- 
portional to muon energy loss on atoms of the mixture 
components. Actually, there is no simple relation be- 
tween the stopping power of a particular type of atom 
and the probability for muon capture by this atom. 
For the binary mixture He/H it is convenient to 
express the atomic capture probalities Wne , Wr by 
per-atom capture ratio A (reduced ratio) defined as 
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where c = CHelCn is the ratio of the helium and hy- 
drogen atomic concentrations. 



probabilities of (7r~)-mesons by hydrogen and helium- 
3 atoms in the H/He mixture, have found that ex- 
pressions ^ fit well the experimental data when the 
slowing-down parameter S is used instead of the cap- 
ture ratio A 
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Here S ~ shc/sh is the ratio of the averaged stopping 
powers, where 



He, H, 



(3) 



are per-atom stopping powers of '^He and H atoms (ex- 
pressed in MeV • cm^/atom), ^ is the mass thickness 
and Ai are the atomic masses. 

Formulae Q can also be used for negative muons cap- 
tured in the H/He mixture (with a proper S*- value) due 
to similarity of fi~ and tt" masses (207TOe and 273me, 
respectively) . 

Recently a series of experiments on the study of fi- 
atomic and /x-molecular processes in a D/'^He mixture 
has been carried out at the Paul Scherrer Institute me- 
son factory |2l I23L |23 | . In the experiments we mea- 
sured the following characteristics: the nuclear fusion 
rate in a charge-asymmetrical muon complex (d/i'^He); 
the probability for transition of the dfj, atom from the 
excited to the ground state (gis); intensities of x-ray 
radiation of /.iHe atoms resulting both from muon cap- 
ture by '^He atoms and from transfer of the muon from 
the d/i atom under its de-excitation to the '^He nucleus. 
For correct interpretation of the above-mentioned ex- 
perimental data it was necessary to have information 
on the probability for direct capture of muons by deu- 
terium and helium atoms in the D/'^He mixture. On 
the one hand, the value of the capture ratio A averaged 
over the data of the papers HIlllllSllllHIll can 
be used; on the other hand, an attempt may be made 
to get experimental information on S = Sjie/ d as a ra- 
tio of stopping powers of helium and deuterium atoms 
in an independent way from an additional experiment. 
The description of the method and analysis of the re- 
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the ratio of stopping powers of helium and deuterium 
per atom is the aim of this work. 



number of electrons detected by the electron counters 
will change when the stopping distribution is shifted. 



II. MEASUREMENT METHOD 

In our experiment we used two types of gas targets; one 
was a D + 5%'^He mixture with the density tpmix = 
0.0585 (hereafter all atom number densities are nor- 
malized to the liquid hydrogen density, Uo = 4.25 • 10^-^ 
cm~'^), the other was a target with pure helium-3. A 
set of different helium densities (pHe was considered. 
The deuterium-helium mixture and the helium targets 
were exposed to the same muon beam in order to keep 
the same initial energy distributions of muons entering 
the targets. The momentum of the muonic beam was 
chosen such as to stop all entering muons inside the 
D/^He target. In the final stage of the slowing-down 
muons are captured by the target atoms, form muonic 
atoms and finally decay via the n~ e~ + + 
reaction. Decay electrons are then markers of the stop- 
ping events. A schematic view of the experimantal 
setup is shown in Fig. ^ (see |2^ and the next 
section for details). 
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FIG. 1; Scheme of the experimental setup. E 
counters. Si - silicon telescopes. 
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Varying the gas density in the pure helium target we 
changed the position of the maximum in the muon 
stopping distribution along the target length. The ex- 
ample of the Monte Carlo simulations presented in Fig. 
121 illustrates such a situation. The stopping code |26l | 
with Ziegler parameterization of the stopping powers 
|27l | was used for obtaining the muon stop distribu- 
tions in the targets. Calculations were performed for 
the 34.0 MeV/c muon beam. Stopping distributions 
were calculated for three targets: the D-|-5%'^IIe mix- 
ture with the density ^Pmix = 0.0585 (thick solid line), 
pure ^He with density (^j^e = 0.0380 (dashed line) and 
pure "^He with the density ifji^ ~ 0.0330 (dotted hue). 
Vertical lines show the position of the central (spher- 
ical) part of the target; this part is directly seen by 
the electron counters. The stopping coordinate a; is a 
distance to the stopping point taken from the entrance 
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FIG. 2: Stopping distributions calculated for three tar- 
gets with different densities ip: D/^He mixture (}Pmix = 
0.0585, thick solid line), ^He((^He = 0.0320, histogram (a)), 
^He((^He = 0.0380, (b)). The beam momentum is 34.0 
MeV/c. The energy of incident muons, is shown in 

the insert. 

The density of the pure helium target can be chosen 
such that the number of muon stops detected via de- 
cay electrons is the same as in the D2 -|- '^He mixture 
(we denote this density by pHe)- In such a case spa- 
tial distributions of muon stops are equivalent. Such 
equivalency of the stopping distributions was verified 
by performing the MC simulations with a set of differ- 
ent ^He densities and examining the differences by x^- 
analysis for both distributions (from D/'^He and '^He 
targets). Figure O shows, as an example, the distribu- 
tions where the minimum of x^/df = 0.92 was achieved 
for pure helium density ip>He ~ 0.0342. Similar equiv- 
alence of the stopping distributions is obtained in the 
plane perpendicular to the beam axis. The details of 
the experimental results discussed in the next section 
also justify the assumption of equivalence of spatial 
distributions of muon stops for a selected He target 
density. 

From the identity of the spatial distributions in both 
targets (considering that the initial muon energy distri- 
butions are also identical) follows equality of the ranges 
of stopped muons for any initial muon energy Ein 
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The above equation can be rewritten in terms of the 
atomic stopping powers s^e, sd 



dE 



dE 



(5) 



3 




2 4 6 8 

stopping coordinate, x (cm) 



FIG. 3: Calculated stopping distribution for the D/^He 
target (histogram, solid line) and the equivalent distribu- 
tion for pure the ^He target with the density (pHe = 0.0342 
(open circles). 



where S is the ratio of stopping powers being the sub- 
ject of the measurement 
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In appendix ^ we argue how a simple and more use- 
ful formula for the mean ratio of the stopping powers 
can be derived from relations Q,© when the behav- 
ior of the individual stopping powers of helium-3 and 
deuterium is taken into account. Such a formula reads 



The experiment included four runs with pure '^He and 
one run with the D/'^He mixture. The experimental 
conditions are given in Tabled 



TABLE I: Conditions of the experiment. Chc is the atomic 
concentration of helium, A'^^ is the number of muons that 
entered the target. 



Run 


Target 


Temp. 


Presssure 












[K] 


[atm] 


[LHD] 


[%] 


[10«] 


1 


3He 


32.9 


6.92 


0.0363 


100 


1.3625 


2 






6.85 


0.0359 




0.7043 


3 






6.78 


0.0355 




0.7507 


4 






6.43 


0.0337 




0.4136 


5 


D/^He 


32.8 


5.11 


0.0585 


4.96 


8.875 



Information on the distribution of muon stops in the 
target volume in experiments with pure '^He and with 
the D/'^Hc mixture can be gained from analysis of 
time distributions of muon decay electrons in '^He and 
D/^He targets 



dN, 
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where E is the average energy of the initial muon en- 
ergy distribution. The above formula gives the recipe 
for measurement of the mean ratio of the helium-3 and 
deuterium atomic stopping powers. 
The specific density (p^e (needed for obtaining S by 
formula Q) is experimentally established by measur- 
ing the yields of electrons from muon decays in the 
D/^He mixture target and in a set of pure '^He targets. 



III. MEASUREMENT AND RESULTS 

The experiment was carried out at muon channel ^E4 
at the Paul Scherrer Institute meson factory. An ex- 
perimental setup (see Fig. developed for studying 
the muon-catalysed nuclear fusion reaction d^^Hc 
a + fj, + p (14.6 MeV) was used to measure Shc/d- 
The body of the cryogenic gas target was made of pure 
Al in the form of a sphere 250 cm'^ in volume. There 
were five kapton windows 55 to 135 ^ thick in the target 
body. The entrance window for the muon beam was 45 
mm in diameter, its kapton was pressed with a stainless 
steel flange with a 1-mm-thick gold ring inserted in 
it. The other four windows were arranged in a circle 
and were intended for detection of charged products 
of fusion reaction in the d^^}lc reaction and muonic 
X-rays from /iHe atoms. Electrons from the decay of 
muons stopped in the target were detected by four pairs 
of plastic scintillator counters (E) installed around the 



dt 



(9) 

where Baz, Bau, Bhc, Fai, Fau and Fd /a He are nor- 
malized amplitudes, B and F are the levels of acciden- 
tal coincidences, and Xai, Aau, ^He, ^o/^He a-re the 
rates of muon disapperance in the target wall material 
and in the target gas. 

Measuring the normalized partial amplitudes B^e and 
Fo/^He and knowing the muon decay electron detec- 
tion efficiency averaged over the muon energy distribu- 
tion, we can determine the number of muon stops in 
^He and the D/'^He mixture. 

Comparing the results of the analysis of the data ob- 
tained under different experimental conditions we in- 
troduced a quantity R for convenience, i? is a ratio 
between the number of electrons from decays of muons 
stopped in ^He (or in D/'^He mixture) and the number 
of incident muons 



R 



1% 



(10) 



The scenario of the experiment was as follows. First, 
the target was filled with a D /^He mixture at the den- 
sity tpmix ~ 0.0585 and the initial muon beam momen- 
tum was varied to find its value corresponding to 
the maximum R. This value of R corresponded to the 
highest density of muon stops in the gas at the cen- 
ter of the target. Then the target was filled with pure 
■^He to the pressure at which the number of stops in 
the target approximately corresponded to the number 
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exact, at which the numbers Ne and iV^ of de- 
tected electrons from decays of muons stopped in the 
gas (in ■^He and D/'^He mixture) were rather close. The 
^He pressure in the filled target was selected using the 
program for calculation of muon energy loss on the pas- 
sage of the muons through various materials [26l |. 
To find out how the ratio R varies with the ^He den- 
sity, four runs were carried out with the helium density 
varying in the interval 1% -7% from run to run (see 
Table D). 

To find the number of muon stops in runs with '^He and 
the D/^He mixture, time distributions of muon decay 
electrons were approximated by expressions © and ^ 
(see ref. 13 for more detailed description). Table ITTl 
presents the numbers of detected muon decay electrons 
as well as the ratios R ~ Ne/N/^i measured in runs 1-5 
with pure helium and D/^He mixture. The results of 
the measurement are also shown in Fig. ^ where four 
experimental points of R for pure helium versus target 
density are plotted. The abscissa of the intersection 
point of i?((/Ji/e)-dependence with the horizontal line 
representing the R{(pmix) value for the D/'^He mixture 
is just the one to determine the value of the equivalent 
helium density (fine- 

In order to evaluate the character of -R-dependence on 
ipHe the auxiliary Monte Carlo calculations were per- 
formed for our experimental conditions. The result of 
the simulation (dash-dotted line in Fig. 0)) shows the 
non-linearity of the R-dependence. 
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FIG. 4: R=Ne/ Nfj, is the ratio between the number of 
decay electrons detected by the electron counters and 
the number of incident muons (black circles) for four 
pure '^He targets used in the experiment (densities: — 
0.0337, 0.0355, 0.0359, 0.0363). The horizontal solid line 
represents the experimental value of R for the D /''He target 
(mixture of deuterium and 5% of "'He, ip = 0.0585). The 
dash-dotted line is the Monte Carlo simulation normalized 
by the factor Rexpifmix) / RMci^Pmix) ■ MC calculations 
were performed using theoretical helium stopping powers 
scaled by the factor of 0.96 in order to obtain better de- 
scription of the experimental points. 



From the analysis of the data presented in Fig. 2|the 
equivalent helium density was found: (p[j^ = 0.0363 ± 



value of the stopping power ratio of helium-3 and deu- 
terium atoms is 



Sne/D = 1.66 ±0.04. 



(11) 



As was mentioned earlier, the fact of the identity of the 
muon stop distributions for the D/^He target at ipmix 
and the He target at (pue is crucial for our analysis. 
There are noteworthy points confirming this fact: 
(i) The numbers of muon stops in the entrance ring of 
Au and the target walls of Al are equal (per incident 
muon) in both cases; 

(m) The ratios between the numbers of stops in the 
target walls and the gas are equal in both cases: 



An additional remark concerning the interpretation of 
our measurement can be made. It is clear that the 
muon slowing-down and subsequent atomic capture are 
closely connected to each another. A discussion pre- 
sented in Appendix ^ shows that such connection is 
so close that allows us to enlarge the interpretation of 
the measurement result from the point of view of the 
atomic muon capture. Concluding, we state that the 
result of the experiment in question is a ratio Sue/D 
of specific stopping powers of helium and deuterium 
atoms while, on the other hand, this quantity is phe- 
nomenologically equal to a reduced ratio of probabil- 
ities, A, for muon capture by helium and deuterium 
atoms. 



APPENDIX A: MEAN RATIO OF STOPPING 
POWERS. 



Using notations 
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we can rewrite the equality of ranges lO as 



v{E) G{E) dE = G=1 
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Two components of the integrated function in Eq. (|A3() 
depend on the muon energy E in a, quite different man- 
ner as is seen in Fig. [3 

The weighting function i'{E) (normalized to unity) is 
strongly energy dependent and decreases roughly lin- 
early with decreasing energy. G{E) (energy dependent 
via S{E)) is, contrary to v{E), approximately constant 
in a wide energy region. In Fig. |B1 the calculated ra- 
tio of the muon ranges in D/'^He and pure helium-3 
targets G = /g i^iE) G{E) dE as a function of the 
initial muon energy, i^m, is presented. For given den- 
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TABLE II: Ratio R measurements for runs 1-5. Nfj, is the number of muons entering the target and is the immber 
of detected electrons from muon stops in the ''He and D/'^He targets. 



Run 


Target 








7? 






[LHD] 


[10«] 


[10«] 


[io-«] 


1 


3He 


0.0363 


1.3625 


0.5302 (14) 


0.3891 (10) 


2 




0.0359 


0.7043 


0.2765 (10) 


0.3926 (14) 


3 




0.0355 


0.7507 


0.2975 (10) 


0.3963 (14) 


4 




0.0337 


0.4136 


0.1657 (8) 


0.4007 (18) 


5 


D/^He 


0.0585 


8.875 


3.4635 (35) 


0.3903 (4) 
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FIG. 5: The functions iy{E) and G(E) (formulae (ETJ and 
IIA2II . respectively). G{E) is calculated for the densities 

IfiHe = 0.0342, iprm^ = 0.0585. 
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FIG. 6: Ratio of muon ranges in D/'^He mixture and pure 
helium-3 targets G = /q i^{E) G{E) dE calculated for dif- 
ferent initial muon energies Ein, for three helium-3 target 
densities: 0.0335, 0.0342, 0.0349. The real energy spectrum 
of muons entering the targets is also shown (top). 



of Ein ■ It is a consequence of the behaviour of the in- 
tegrated function (or, in other words, due to similar 
dependence of deuterium and heUum stopping powers 
on the muon energy). In the energy interval 1.8 - 2.3 
MeV (50 % of beam muons belong to this interval) the 
relative change of G is 0.4 %, and for interval 1.3 - 2.7 
MeV (90 % of muons) the respective change is 1 %. 
The quantity G in Eq. IjASp represents the ratio of the 
ranges of muons with the initial energy Ein in D/'^Hc 



is fulfilled for a given energy Ein for especially 
chosen (pHe (as is seen in Fig. 0. For another energy 
the other density (pHe should be, in principle, adjusted. 
But as can be seen from Fig. |3 such uncertainty in (pue 
is very small (less than 1 % in the range of our muon 
energy spectrum) and can be neglected. 
In view of the above considerations it is reasonable to 
use an approximation 



G « G{S{E)) w G{S{E)), 



(A4) 



where S{E) is the ratio of the atomic stopping powers 
taken for the average energy of the initial muon spec- 
trum ^2 MeV. Then from cqs. l|X3)l and (|X4jl follows 
the equality 



1, 



(A5) 



and finally formula 0. 

APPENDIX B: ATOMIC CAPTURE 

Slowed-down muons are captured by the target atoms 
and finally end their life decaying in the atomic orbit 
0,123. ^ 

An alternative phenomenologocal interpretation of our 
measurement is the following. Each muon stop in the 
target gas is followed by an atomic capture cither on a 
deuterium atom or on a helium atom. We identify the 
stopping event with an atomic capture event because 
we detect the decay electrons appearing after the muon 
is captured (electron decay in flight can be neglected 
as the slowing down time is about 1 ns, the probability 
of such decay is ^ 10""'). The number of decay elec- 
trons, which appear as the consequence of muon stop 
at a given point of the target, depends on the density 
of atoms and on the probability of the atomic muon 
capture. 

Let iVp muons enter the target with a given initial en- 
ergy distribution. Through a volume located at 
any given point of the target there pass ATV^ muons. 
A part of them, AA^^, having a sufficiently low en- 
ergy due to slowing-down, arc captured in AT^ and 
produce decay electrons. These muons which arc not 
captured in AV continue to slow down and arc stopped 
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electrons (i.e. the number of stopped and then cap- 
tured muons) in AV^ is proportional to the number of 
entering muons AA^^ and to the muon capture proba- 
bility described by the capture cross section ctc. The 
mean muon capture probability in AV is ScAx, where 
Ec = n{He)aHe + n{D)aD is the macroscopic capture 
cross section, n are the densities of helium and deu- 
terium atoms, WHe and are the muon capture cross 
sections for helium and deuterium atoms, respectively 
(averaged over the muon energy), and Ax is thickness 
of the volume Ay along the direction of the muon 
beam. Under the assumption that the muon atomic 
capture cross sections aneiE), aoiE) arc the same in 
both cases (monoatomic or mixture target), and do not 
depend on the He concentration, the respective num- 
bers of the captured muons can be written as 



as a general relation 



AiV,^ 



stop 



= AN^UoipHef^HeAx, 



(Bl) 



ANZ; = AN^Uof^mACD-^D + CHeaHe)Ax, (B2) 

for a pure helium-3 target and for a target with a 
D-|-'^He mixture, respectively. 

For a chosen density of the pure helium target, pHe, 
the spatial stopping distribution is the same as for the 
mixture target (such a situation is illustrated in Fig. 
13). Then, for the same volumes AV in both targets, 
the number of stops in helium the target with density 
PHe is equal to the number of stops in the mixture with 
the density pm,ix 



AN"t^ = ANT"" 

'—^^^ stop stop- 



(B3) 



Since the above local equality is valid for any point of 
the target we can rewrite it (using cqs. (|B1|) and (|B2p 'l 



where 



a{He) 



(B4) 



(B5) 



is the ratio of the probabilities of muon capture by 
helium-3 and deuterium atoms. Rewriting Eq. ljB4p in 
the form 



A = 



(B6) 



one obtains the formula for the quantity A, which is 
equivalent to Eq. 

From the above considerations and taking, in particu- 
lar, into account the identity of the relations ||7J) and 
(jB6ll . one can state that the value of the atomic stop- 
ping power ratio S, which we measure in our experi- 
ment, is also (under indicated restrictions) a value of 
the reduced probability ratio, A^ for the muon capture 
by helium-3 and deuterium atoms. 
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